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Abstract 

We calculate the distributions in recoil momenta and their energy distributions for the high 
energy non-relativistic double photoionization of helium caused by the quasifree mechanism of the 
process. The distributions obtain local maxima at small values of the recoil momenta. This is in 
agreement with the earlier predictions and with recently obtained experimental data. We obtained 
also the angular correlations, which reach the largest value in the "back-to-back" configuration 
of the photoelectrons. Our analysis is valid in all high energy nonrelativistic region. Particular 
equations are true for the case when the wavelength of the photon exceeds strongly the size of the 
atom. We present numerical results for the photon energies in the region of 1 keV, employed in the 
recent experiments. 

PACS numbers: 32.80.Fb, 34.80.Dp, 31.15.V- 
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I. INTRODUCTION 



?he recent measurement of the yield of double charged ions in photoionization of helium 
3| confirmed the exitance of the quasifree mechanism (QFM) of the double photoioniza- 



tion which was predicted in [4j. The differential 
were calculated earlier in a number of papers 



cross sections of the double photoionization 
4 ] — (9 1 , where the authors studied the distri- 



butions in characteristics of each photoelectron. In the pioneering experiments the 
distribution in momentum transferred to the nucleus q (recoil momentum) was measured. 
Thus the problem of calculation of such distributions as da 2+ /dq 2 de with e the energy of 
one of the photoelectrons and da 2+ /dq 2 became actual. 

In the present paper we calculate these differential cross sections for the double pho- 
toionization of helium at high values of the photon energies, corresponding, however, to 
nonrelativistic energies of the photoelectrons. We present the results for the differential 
cross sections da 2+ / dq 2 de We trace the dependence of these characteristics on the photon 
energy uj. 

Recall that the experiments in which the photons carried the energies 450 eV, 800 

eV and 900 eV demonstrated that the distribution of outgoing electrons obtains a surplus 
at small q of about 2 a.u. The kinematics of these experiments enables to separate the 
non-dipole contributions at small values of q. Thus the observed surplus is entirely due to 
the non-dipole terms. 

By that time only two mechanisms of the process were known. In both of them the 
electron which interacted with the photon directly obtained almost all the incoming photon 
energy uj. In the first one, known as the shake-off the secondary electron is pushed to 
continuum by to the sudden change of the effective field. In the second, called the knock-out 
mechanism, the photoelectron inelastically collides with the bound one, sharing the photon 
energy. Both mechanisms contain the single photoionization as the first step. This process 
can not take place on a free electron. Thus the momentum q which is transferred in this 
step to the nucleus exceeds strongly the averaged momentum of the bound electron /1, in 
the case of the high energy photon 

w>/, (1) 

(J is the single-particle ionization potential). Of course, there is a configuration in which 
the second electron transfers momentum qi 3> [i such as |q + qi| ~ fi- However since each 
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act of transferring a large momentum g>/j leads to an additional small factor (lOj, [ll|] its 
probability is very small. The distribution da 2+ /dq 2 de provided by these two mechanisms 
peaks at q ~ (2mw) 1 / 2 > /i (m is the electron mass), becoming very small at q ~ //. This 
remains true beyond the dipole approximation. 

In contrast to a single electron case, the two electrons can absorb a photon without 
participation of the nucleus. In the free process q = 0. In the QFM small momentum 
q ~ /i is transferred to the nucleus, i.e. q is much smaller than the momenta of the outgoing 
electrons. The distributions da 2+ ldq 2 de and da 2+ Idq 2 the have local maxima at small q of 

n n 

the order of rj. That's what was detected in [II] -[3J. 

Momentum q transferred to the nucleus can be written as 

q = k-px-p 2 , (2) 

where pi j2 are momenta of the outgoing electrons, while k is that of the photon. The recoil 
momentum q can become small only if the large momenta of the outgoing photoelectron with 
Pi — I Pi I yU compensate each other to large extent (k = |k| is always much smaller than 
Pi while we consider the photon energies, corresponding to nonrelativistic photoelectrons) . 



Hence the values of pi should be close, i.e. p± ~ P2 ~ yrnE with E the sum of the energies 
of the photoelectrons. Thus in QFM the bound electrons exchange by small momentum 
q ~ fi with the nucleus and by large momentum of the order pi ^> /i between themselves. 

We calculate the amplitude of the QFM in the lowest order of expansion in powers of q/pi- 
This corresponds to expansion of the bound state wave function in the lowest order in powers 
of r 12 /rj, with standing for the distance between the electron and the nucleus, while r 12 
is the interelectron distance. We consider the high energy photons, corresponding however 
to nonrelativistic energies of the outgoing electrons. Thus we assume that u <C m. Having 
in mind future extension of the analysis to the relativistic case we employ the relativistic 
units h — c— \. 

Since q ~ fi the higher terms of expansion in powers of qjpi are of the same order as those 
coming from the interactions between the photoelectrons and the nucleus. However, they 
are of quite different physical origin. Thus we include interactions of the nucleus with the 
electrons exactly, describing the latter by the nonrelativistic functions of the Coulomb field. 
Interaction between the outgoing electrons is proportional to the square of its Sommerfeld 
parameter £ ee = a/v, where v is their relative velocity. For the energies of the order 1 keV, 
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employed in the experiments 

Q-S 

this interaction provides a correction of the order of 2 

% and can be neglected. 

Direct relation of the QFM to the behavior of the bound state wave function ^(ri, r 2 , r 12 ) 
at small distance r 12 was demonstrated in [loj]. It was shown that the QFM amplitude 
contains the factor <9\I//<9r 12 at r 12 = 0, which is connected to the function \I/(ri, ?"!, r 12 = 0) 



by the Kato cusp condition The latter appears to be very important for calculation 
of the QFM amplitude p|. We employ very precise wave function [3] which satisfy also 
the Kato cusp conditions. We use an analytical function which approximate these wave 



13| 



functions at the electron coalescence line r 12 = very accurately 

We include only the quadrupole part of the electron-photon interaction. Note that the 
amplitude contains also the dipole terms proportional to the product (e-q) with e the vector 
of the photon polarization. However, at least in the leading approximation it is canceled 
by the contribution in which the electrons exchange by large momentum in the final state 
7j. Anyway, in the experiments ;he observations were carried out in the plane where 

(e ■ q) = 0. Thus we can focus on the quadrupole contribution. 

Besides the conservation of the linear momentum expressed by Eq.((2]) we write the energy 
conservation condition 

u-I = E; E = e l + e 2l (3) 

where £j = pl/2m {i = 1, 2) are the photoelectron energies. Note that in our system of units 
cu = |k|. 

To simplify the calculations we restrict ourselves to the case when the photon wave length 
is much larger than the size of the bound state, i.e. 

u < /J,. (4) 

For the atom of helium this means that w C 6 keV. Under this condition Eq.(T5]) can be 
written as 

q=-Pi-P2, (5) 

in the lowest order of expansion in powers of k. 

Momentum q can become as small as n only if momenta of the outgoing electrons almost 
compensate each other, i.e. |pi + P2I ~ H <C Pi,2- Hence, the photoelectrons are emitted 
mostly "back-to-back, with t = (pi • P2VP1P2 close to —1, while the values pi ~ p 2 , i.e. 
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\pi — P2I *C pi,2- Thus the relative difference of the energies of the outgoing electrons 

(6 ) 

should be small. Since q > \pi — Pz\ we find 

This equation is presented in the lowest order in /3. 

Besides the distributions da 2+ / dq 2 de and da 2+ /dq 2 we calculate the differential cross 
sections da 2+ /dtde and da 2+ /dt. We present the numerical data for the photons carrying 



fl-S. 



the energy of about 1 keV, employed in the experiments 

Note that this approach was used in [ijj] for calculation of the distributions da 2+ /dtde 
and da 2+ /de at the point of the peak t = — 1. In other words in [14|] the height of the peak 
of this distributions was found. In the present paper we calculate the shape of the peaks. 



II. GENERAL EQUATIONS 

The differential cross section of the double photoionization can be written as 

da 2+ = ^-\F(k, Pl ,p 2 )\ 2 dT. (8) 

Here -F(k, p 1; p 2 ) is the amplitude of the process. Averaging over polarizations of the photon 
is assumed. The last factor is the phase volume 

dV = 2n5(u-I-e 1 -e 2 ) {27r)3(27rr (9) 

Employing Eq.([5]) we can present 

dT = 5{u - 1 - 2ei - ^ - i-)*tp*p$-. (10) 
v m 2m' 4vr (2tt) 3 v ; 

with z the direction of momentum p 1 . Using 5-function for integration over q z we can write 

2+ . ( 2 6(q/p-P) Em 2 d/3 dQ dq 2 
da =|F( k , Pl) p a )|_ (11) 

with Q the solid angle of the photoelectron with momentum pi, p = [mE) 1 ^ 2 . 
The amplitude of the process can be written as 

F(k, Pl , p 2 ) = 2)| 7l + 72 |^(1, 2)), (12) 
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with the numbers 1 and 2 denote the variables corresponding to two electrons, V&jj are the 
wave functions of the initial and final states, 

7 = (47ra) 1 /V( k - r ) ~' (e ' V) , (13) 

m 

where e is the vector of polarization of the photon, (e-k) = 0. Recall that we shall pick only 
the quadrupole terms of interaction between the photon end electron. For further evaluation 
we denote 

F(k, Pl , p 2 ) = (4™) 1 / 2 M(k, Pl , p 2 ). (14) 
As we said earlier, we describe the final state by the function 

*/(!-!, r 2 ) = -J=^i(r 1 )^ 2 (r 2 )+^2(ri)^i(r 2 )), (15) 

where ipi are the single-particle nonrelativistic Coulomb field function with asymptotic mo- 
menta pi. We shall need the functions 

^ i (r) = e-^- r )X( Pi ,6,r); » = 1, 2 (16) 

Here 

X( Pi ,&,r) = iV^OiFx^i, l,i Pi r + i( Pi • r)), (17) 

£i = — ; i] = maZ. (19) 

Pi 

Here Z is the charge of the nucleus. Note that in the hydrogenlike approximation rj is the 
averaged momentum of the electron in the Is state. 
Thus 

M(k, Pl , p 2 ) = V2(A(k, Pl , p 2 ) + ( Pl <— > p 2 )) ; (20) 

A(k, Pl ,p 2 ) = J rf 3 r 1 rf 3 r 2 e i ( k ri )- i ( pi ri )-' i ( p2 r2 )X 1 (r 1 )X 2 (r 2 )(e- V ri )^( ri ,r 2 ) = A. 
Here we denoted 

X 4 (r)=X( Pi ,&,r). (21) 

Introduce 

R=^±^; p = r 2 -r 1 (22) 



while 



with 



Presenting 

ri = R - p/2; r 2 = R + p/2; V n = ^V fi - V p , (23) 

and 

^(n, r a ) = *(R,p) (24) 

we obtain 

A = A! + A 2 , (25) 

with 

A 1 = - [ rf 3 i?d 3 pe-^ +i ( q R )X 1 (R-^)X 2 (R+^)(e-V p )v&(R,p), (26) 



while 



A 2 = -^| ^V-* (a - p ^^ (27) 



where 

Pi - p 2 + k 



a = |a| (28) 



2 

Since a 3> g, the integrals on the right hand sides of Eqs.()26l). (12T1) are saturated by 
i? ~ l/g> p. Thus we can put p = in the functions X;. This provides 

A x = ^-^ /" d 3 Re i ^ R) X 1 {K)X 2 {K) J d 3 pe^ (a ' P) ^(R, p), (29) 

Now we expand the wave function 

*(R, p) = r, p) = *(f2, 0, 0) + tK(R, r, 0) + pV' p (R, 0, p) + 0(p 2 ). (30) 

Here r = (R • p), the derivatives are taken at r = p = 0. 

We can calculate the integral over p multiplying the integrand by e~ vp and putting v — 
in the final step. 

/ ,V^P)*(R, p) = / ^(a-P)-^ (R) p=0) = / ^(a-P)-^^ 0) , = 0) = 

(31) 

8^(^,0, p = 0) 
a 4 

The derivative \l/ p at p = is related to the wave function by the Kato cusp condition 
13| 
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lim p -> r y' p (R, p) = ^(R, 0), (32) 

where tq = 1 /ma is the Bohr radius. It is identical to similar relation for the wave function 
presented in variables ri,r 2 ,p. Introducing 

$(i2) = *(R,p = 0), (33) 

we can write 

^^(e-a)^), (34) 



a 4 



with 



with the functions -Xj(R) defined by Eq. (|2ip . 

Combining Eqs. (|T4"t |20| |25| [28]) we find for the quadrupole terms of the amplitude 



S^q) = J d 3 i?e i ( q R )X 1 (R)X 2 (R)$( J R), (35) 



F(k, Pl , p 2 ) = (4vr a ) 3 / 2 4v / 2 (e " n)( 4 k ' n) g 1 (g). (36) 



P 4 



After averaging over the photon polarization and integration over the angles Eq. fllll) takes 
the form 

7^ a Tz\ 3 W)\ ■ ( 37 ) 



dq 2 d/3 15 E 4 

In order to calculate Si (q) we employ the presentation of the function 

$(#) = cie~ Al/? + c 2 e~ A2i \ (38) 
with numerical values of the parameters 

ci = 0.380C 3 , c 2 = 0.990C 3 , Ai = 5.54C, A 2 = 3.41C, C = ma 



found in [13]. 



The further calculations are described in Appendix. We obtain for the function S\(q) 
defined by Eq. fl35|) 

|Si(9)| 2 = l$> J ( A *)| 2 ; * = 1 > 2 - (39) 



with defined by Eq. CATTT) . 
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III. RESULTS 



Now we present the results of computations. The cross section d 2 a / 'dq 2 df3 determined by 
Eg. ( 157)) is presented in a three-dimensional Fig.l for u = 800 eV. As expected, it obtains 
the largest values at small/3 1 and in the region of small q ~ la.u. in agreement with 
the experimental results {3]. This distribution at /3 = 0, corresponding to the center of 
the spectrum is shown in Fig. 2 for u = 800eV and u = 1 keV. Since the effects of finite 
(3 manifest themselves in the terms of the order /3 2 , there is no noticeable difference from 
similar figures for /3 7^ in the QFM region due to Eg.(J7J). 

It is instructive also to view the energy distribution of the angular correlation 

d2 ° o d2 ° + (Pi ■ P2) (An s 
dtdp dq z dp V\Vi 

It is shown for u = 800 eV in the three-dimensional Fig.4. As expected, the largest values 
are reached at (3 1 and t close to —1, corresponding to the electrons ejected in the opposite 
directions ("back-to back"). For /3 = this differential cross section is shown in Fig.4 for 
u = 800eV and u = 1 keV. 

We calculate also the distribution in recoil momentum 

da_ 1 f q/p , R J?o_ 

dqi~2j P dqH^ [ ' 

and the angular correlation 

da 1 f 1 ,„ d 2 a 



dt 2J dP dtd^ (42) 
They are presented in Fig. 5 and Fig. 6 correspondingly. As expected, the distribution da/dq 2 
has a local maximum at q about 1 a.u. At q = this distribution turns to zero just because 
the interval of integration over /3 vanishes. The angular correlation da /dt has a sharp peak 
at t = — 1, in agreement with the previous analysis. 



IV. SUMMARY 



We calculated the distributions in recoil momenta q and their energy distribution for 
the high energy nonrelativistic double photoionization of helium caused by the guasifree 
mechanism (QFM)[4|. They are closely related to the distributions in the angle between 
momenta of the outgoing electrons (angular correlations). As expected, the distributions in 
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recoil momenta obtain local maxima at small q of the order 1 — 2 a.u., in agreement with the 



results of 
results in 



;he pioneering experiments 



. Unfortunately, the way of presentation of the 
l[-|3j does not permit to compare the quantitative results. The corresponding an- 
gular distributions obtain maxima when the photoelectrons move in the opposite directions 
("back-to-back" scattering). The qualitative picture is the same for heavier atoms. 

The QFM is caused by the initial state interactions and, contrary to a misleading state- 
ment in [3|, is contained in the standard Feynman diagrams for the amplitude |4|. Since 
the QFM is at work at small separation between the bound electrons r±2, we described the 
initial state by a very precise wave function H, e,np l0yi n g its ana lyt ica, approbation 
at small values of r 12 13]. We neglected the electron interactions in the final state. The 
numerical results for the photon energies in the keV region are shown in Figs. 1-6. This 
energy region attracts attention nowadays in connection with the laser experiments. Also, 
the experiments where carried out at these energies. The approach can be applied for 

the double photoionization of heavier atoms. 

One can obtain more precise results by direct employing of precise wave functions, i.e. 
those found in However, such approaches do not allow to analyze the mechanisms of 
the process. On the other hand very "accurate" wave functions (i.e. those which reproduce 
the value of the binding energy very accurately) may have a wrong behavior at ri2 — > 0. 
It was demonstrated in 



[l6| that a number of publications on the subject employing 
such functions contain erroneous results. More examples are given in That is why we 



consider our results as a necessary step in investigation of the process. 

There is a number of possibilities to carry out the experimental investigation of other 
phenomena connected with the QFM. Outside the plane (e • q) = interference between 
the dipole and quadrupole terms should manifest itself in the angular distributions. Also, it 
would be interesting to trace the ui dependence of the shape of the energy distribution. Its 
theoretical analysis was presented in 171 ]. 

The work was supported by the MNTI-RFBR grant 11-02-92484. One of us (EGD) 
thanks for hospitality during the visit to the Hebrew University of Jerusalem. 
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Appendix A. 



Thus we calculate the integral 

S 1 {q) = J d 3 Re l ^X 1 (R)X 2 (R)e- XR = N(^)N(^)I(\), (A.l) 

where 

/(A) = J(A) = J d 3 Re^F 1 (R)F 2 (R)^, (A.2) 

with 

F l= xFiti&MPi-R + iCPi-R.)), (A.3) 



The integral J(A) was calculated in [15[ as 



J(A) = £^(2£)*(i±«.)-* 2 F l( l- i6 ,i6,l, 9 ), (A.4) 



with 

g 2 + A 2 

a c = — g — ' ^ = (P2-q)-«Ap 2 ; 7c= -(Pi-q)+iApi-a c ; 5 C = p x p 2 - (pi • p 2 ) - Pc, 

(A.5) 

_ aJc ~ Prfc 
a c (7c + S c ) ' 
We write Eq( 1A.4l) in a more symmetric form 



m = ^^fv * f 2Fl (i^i,hi (a.6) 



(A.7) 



with 

/3 C 7 C - a c <5 c 



7 C (a c + p c ) ' 
Thus 

J(A) =4vrA(A) 2 F 1 (^i,^2,l,MA)), (A.8) 

with 

A(A)=(y + A-J (pi+p 2 + iAj (pa-pi-iA) (pi-p 2 -iAj . (A.9) 

Employing 

2^(^1,^2,1^) = -Ci6 2 *i(i£i + M6 + i,2,/0, (A.io) 



we find 



J ( A ) = (g 2 + A 2 )2 et(6+6) ( A ) T ( A ) e ' V2(€l+fe) - ( A - U ) 
11 



Here 



9(A) 



q 2 + A 2 
s{X)u{\) 



>(A) = v / (pi+p 2 ) 2 + A 2 , u(A) = v / (Pi-P 2 ) 2 + A 2 



while 



T(A) = (1 - !iil±M(i + /i(A)) 2 Fi(z^i, i£ 2 , 1, fc(A))- 
- MA)) 2 Fi(*6 + 1, «6 + 1, 2, fc(A)), 



(A.12) 
(A.13) 



with 



h(\) = 1 



g 2 + A 2 
m 2 (A) 



(A.14) 
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FIG. 2: Distribution da 2+ /dq 2 d/3 in KT 10 ^, r = I J ma for (3 = 0. Solid line is for oj = 800 eV, 
dashed line is for u = 1 keV. The recoil momentum q is in atomic units. 
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FIG. 4: Distribution d<r 2+ / dtd/3 in barns for uj = 800 eV. 
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FIG. 5: Distribution da 2+ / dtd/3 in barns for /3 = 0. Notations are the same as in Fig. 2. 
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FIG. 6: Distribution da 2+ jdt in barns for /3 = 0. Notations are the same as in Fi 
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